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Abstract 

Cosmological models with an SU(2) Yang-Mills field are studied. For a specific model 
with a minimally coupled Yang-Mills Lagrangian, which includes an arbitrary function of 
the second-order term and a fourth-order term, a corresponding reconstruction program is 
proposed. It is shown that the model with minimal coupling has no de Sitter solutions, for 
any nontrivial function of the second-order term. To get de Sitter solutions, a gravitational 
model with nonminimally coupled Yang-Mills fields is then investigated. It is shown that 
the model with non-minimal coupling has in fact a de Sitter solution, even in absence of 
the cosmological constant term. 

1 Introduction 

One of the most important recent results in observational cosmology has been the discovery that 
the Universe is speeding up, rather than slowing down. The combined analysis of data coming 
from Supernovae la (SNe la) pQ, from the cosmic microwave background (CMB) radiation [21 [3], 
large scale structure (LSS) pE], baryon acoustic oscillations (BAO) [5], and weak lensing [6], 
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gives strong evidence for the accelerated cosmic expansion. This acceleration suggests that 
the present day Universe is dominated by a smoothly distributed slowly varying cosmic fluid 
with negative pressure, the so-called dark energy. Contemporary cosmological observation data 
strongly support the conclusion that the dark energy equation of state (EoS) parameter is very 
close to minus one. These observations also confirm that the Universe is, to high approximation, 
isotropic at very large scales. 

Various scenarios for the late-time accelerated expansion of the Universe have been proposed 
(for reviews, see [3, El El HQl [XXI [12] ) . The simplest of them is general relativity with a cosmological 
constant (for a review see, e.g., [13J). The standard way to obtain an evolving EoS parameter is to 
add scalar fields into a cosmological model. The evolution of the Universe is quite well described 
by cosmological models with scalar fields, in particular, by quintom cosmological models, which 
are two-field models with a phantom scalar field and an ordinary scalar one. Quintom models are 
being actively studied at present time [HI [151 EE EL71 EH] ( see a l so P2] f° r reviews). Some other 
popular models involve modifications of gravity, as for instance F(R) gravity, with F(R) an (in 
principle) arbitrary function of the scalar curvature R (for recent reviews see, e.g., [TTl [T2"| I2"U]). 

A possible source of both inflation and the late-time acceleration of the Universe is a nonmin- 
imal coupling between the scalar curvature and the matter Lagrangian [211 [22] (see also [23]). 
Such a coupling may be applied for the realization of a dynamical cancelation of the cosmological 
constant [21]. Criteria for the viability of such theories have been discussed in [25[ [261 127] . A 
coupling between a function of the scalar curvature and the kinetic term of a massless scalar 
field was considered |28|. 



It is known that the coupling between the scalar curvature and the Lagrangian of the elec- 
tromagnetic field arises in curved space-time due to one-loop vacuum-polarization effects in 
Quantum Electrodynamics (QED) [29]. A unified inflation and late-time acceleration model, 
with the electromagnetic field coupling to a function of the scalar curvature, has been proposed 
in [30] by using the analyzing procedure in the electromagnetic field [31] • Inflation driven by 
a vector field has been discussed in [32]. A related scenario is known as vector inflation [33J. 
The idea of vector inflation is to use a vector field as the inflaton, the field which gets non-zero 
background value during inflation and drives the inflationary dynamics. 

Non-Abelian gauge fields are widely used in particle physics and are being actively studied 
in cosmology ^31 ISH I3H1 1101 SD S21 S3]- Note that string compactifications may 

naturally lead to an effective-scalar- Yang-Mills-Einstein theory (plus higher-order corrections). 
Inflationary cosmology and the late-time accelerated expansion of the Universe in a non-minimal, 
non-Abelian gauge theory (the Yang-Mills (YM) theory), in which a non-Abelian gauge (YM) 
field plays a significant role, has been considered in [39], where the authors show that the 
appearance of such non-minimal terms in the early Universe can be compatible with current 
formulations of the YM theory coming from a specific choice for the non-minimal function. Also 
in [39] the cosmological reconstruction of the YM theory has been discussed and a correspond- 
ing algorithm has been proposed. A new variant of the accelerating cosmology reconstruction 
program has been developed in [ID]. 

A remarkable fact is that the SU(2) Yang-Mills field admits an isotropic and homogeneous 
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parametrization by a single scalar function. This parametrization, which is useful for the re- 
construction program [39l |4"0] , has been employed to get an inflationary scenario [¥!] . As noted 
in [H] , the standard YM term, minimally coupling with gravity, does not lead to inflation, and 
thus one should add new terms in order to get a convenient inflationary scenario. In [JTj 



an in- 



flationary scenario, in which slow-roll inflation is driven by a non-Abelian gauge field minimally 
coupled to gravity has been proposed. To achieve this, the authors add a fourth-degree term to 
the YM Lagrangian. 

In this paper, we generalize the reconstruction program in the case of a YM Lagrangian with 
a fourth-degree term (see action ([1]) below). We then get the equation which connects the YM 
field with the Hubble parameter. This equation does not depend on the form of the arbitrary 
function J 7 . Thus, assuming an explicit form for the Hubble parameter as a function of time, we 
can get the YM field and, after this, construct the function J 7 , which corresponds to the given 
solution. One of the main results obtained here is that the model with minimal coupling has no 
de Sitter solution, whatever the nontrivial function J 7 be. 

As is well known, de Sitter solutions play a very important role in cosmological models, 
because both inflation and the late-time Universe acceleration can be described as a de Sitter 
solution with perturbations. In order to obtain these solutions, we then consider a gravitational 
model with nonminimally coupled YM fields. We succeed in deriving a first-order ordinary 
differential equation for the YM field which corresponds to the de Sitter solution. This equation 
includes an arbitrary parameter. Depending on the value of this parameter, it can be solved in 
quadratures or else numerically. 



2 Gravitational models with the Yang— Mills fields 



Let us consider a minimal gravitational coupling of the SU(2) Yang-Mills field in the general 
theory of relativity, which is described by the following action: 



S 



d x \/—g 



Ml 



R + F{Z) + 



—A 

384 



r a/3 r Xa) 



A 



(1) 



where g is the determinant of the metric tensor g^, R is the scalar curvature, Mp = Mpi/\/8n, 
the Planck mass M P1 = 1.2 x 10 19 GeV, and k is a constant. The SU(2) YM field has the 
internal symmetry index a, the field strength tensor being 



d a A% - d p A a a + f abc A b a A%. 



(2) 



The function F{Z) is an arbitrary function of Z = F^ u F a ^ v (summation in terms of the index a is 
understood), while the numbers f abc are structure constants and thus completely antisymmetric. 
For the SU (2) group, 

f abc = - ~g[abc], (3) 
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where g is a constant and [abc] the Levi-Civita antisymmetric symbol (we use this notation 
instead of e abc because we reserve the last one for the Levi-Civita antisymmetric tensor). Roman 
indices, a, b, c, will run over 1, 2, 3, and the Levi-Civita tensor is given by 



-99 Pia 9 P ^9 P3X 9 P4CT [piP2P3PA], [0123] = 1. 



(4) 



Models of this kind, in the case R — 0, have been considered in jlQ]. The case TiZ^) = Z has 
been analysed in [UJ where an inflationary scenario, in which slow-roll inflation is driven by a 
non-Abelian gauge field minimally coupled to gravity, has been proposed. 
The equation of motion for the field A a ^ is 



j'{Z) F a ^ + Je^F* 
192 



~g^g-[abc]A b A T' (Z) F 



(5) 



where J = e a ^ Xa F h a gF\ a . Variation of ^ with respect to g^ v yields the field equations: 



M 2 P 



- \g,uR\ - \g^T (z) + it' (z) f; p fz»- 



384 [2 

In the spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) universe 

ds 2 = — dt 2 + a 2 (t) (dx\ + dx\ + dxl) , 
the following ansatz for the SU(2) field: 



A b 



{o,<p(t)s b ). 



(6) 



(7) 



(8) 



is very useful [39l HQl HT]. The ansatz identifies the combination of the Yang-Mills fields for 
which the rotation symmetry violation is compensated by the gauge transformations. Thus, we 
get the rotationally invariant energy-momentum tensor of the Yang-Mills fields, namely 



F b 



and 



F b 



F b - 

1-3 



[bij], F b0j 



-h-8 bj , F blj 
a 2 



Z = F b „F b ^ = 6 



~9 2 ^ 



er a" 

where differentiation with respect to time t is denoted by a dot. 



■[bij] (9) 



(10) 
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Use of the ansatz ([8]) allows to obtain the YM energy-momentum tensor, having the same 
form as an ideal isotropic fluid with the energy density p and the pressure P, in other words: 



T 



(YM)/i 



diag(-p, P,P,P). 



Taking into account (J4]) and ([7]), the equation of motion ([2]) can be written as follows: 



-9 [abc] A b u Lf^F' (Z) F c ^ + J^JfrvaflF^ = 0. 
It is convenient to write the Friedmann equations in terms of ip = <p/a. Using 



a ( ip + Hip 



a [ip + 2Hip + ip(H + H 



where H = a/a is the Hubble parameter, we get the equations which follow. 
The (0, 0) component of (T5]) reduces to: 

3 -^H 2 + X -T (Z) + 6T ' (Z) + H^) 2 - 5«SV 4 (^ + Hip) 2 -\a = o, 
where ip = <p/a. Note that 

Z = 6 (ftp 4 -(tp + Hip) 2 ) , Z = 12 (2g 2 ip 3 ijj - (ip + Hip){ip + + Hip)) , 



(12) 



(13) 



J = 2A~gip 2 [ip + Hip\. 



The (i, i) components of (T5]) yield: 



-M 2 P 



2H + 3H 2 



-T(Z) + 2T'(Z) 



ip + Hip) -2g 2 ip 4 



(14) 



h~g^ + HlP) 2 -h = 0. 



By subtracting Eq. (fl3|) from Eq. (fT5|) . we get 



A/f 2 / / \ 2 

-fH = 2T'{Z)[U + Hip] +~g 2 ip 4 



(15) 



(16) 



From this equation, it follows that the model considered does not have nontrivial de Sitter 
solutions (H is a constant). Such solutions can exist only if either T{Z) is a constant, or the 
function ip(t) = 0. In the next section we show that nontrivial de Sitter solutions do exist in a 
model which has a non-minimal coupling. 

Rewriting the last equation in the following form 



T'(Z) 



M 2 



H[[iP + HiP) + g 2 iP 4 



(17) 



introducing (ITT)) into Eq. (TT3T) and differentiating, we get: 



H 



2[[i> + Hip) + g 2 ip A 



2# W - [ip + Hip) [ip + Hip + Hip)) + 



H 



ip + HipJ + g 2 ip A 
2 • 1 



+ 



hU + Hip\ U + Hip + Hip 



ip + Hip) + g 2 ip A 



- -^f \f + H?P) V + (V> + #V>J (V> + #V> + ^jj = 0. 

If we assume, or rather know, the specific form of the Hubble function H(t), then Eq. ( |T8l) 
constitutes a differential equation for ip(t) and, once we determine this function, the correspond- 
ing Yang-Mills theory can be found (i.e., the function iF (Z)) which reproduces the cosmology 
given by H(t) in the frame of the spatially flat FLRW universe. From (ITTj) we can find J 7 (Z) 
up to an integration constant. This constant can be determined from (fl5|) and corresponds to 
the cosmological constant. 



3 Non- minimal gravitational coupling with the Yang 
Mills field 



3.1 Action and equations 

In this section we will consider a non-minimal gravitational coupling of the SU(2) YM field in 
general relativity, which is described by the action: 



S, 



GR 



d A x\f^g 



S R + £ym — A 



Cym = --(l + f(R))Z, 



(19) 
(20) 



where f{R) is an arbitrary, thrice differentiable function of R. 

The field equations can be derived by taking variations of the action in Eq. (TTD"]) with respect 
to the metric g^ and the SU{2) YM field A^, as follows: 



R^iu c^p^R ~ ftnv ^ ^9 Liu) ? 



(21) 
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with 



+ - {f{R)FR^ + 9ilv U [f{R)T\ - V„0„ [f'(R)T]} 



(22) 



where the prime denotes derivative with respect to R, V M is the covariant derivative operator 
associated with g^, and □ = g^V '^d v is the covariant d'Alembertian for the scalar field. 
It is convenient to write down the trace equation 



R=~ ^f v (^ M) - 9,A) =-^2 {fiR)™ + 3D \f{R)F\ - 8A} 



(23) 



We will show that the trace equation is useful in order to find the de Sitter solutions. 



3.2 The Friedmann— Lemaitre— Robertson— Walker metric and equa- 
tions of motion 

Using the ansatz (jSJ), we get the following equations in the FLRW metric (see the Appendix, for 
details): 



3H' 



M 



1 / . » A 3 
-(A + p) = — + — 



M 2 2M 2 u 



l + f(R)) [g'r+^ + H^Y - 



6 (H + H 2 ) f'(R) g 2 ^ - {i> + Hip) + 6Hd 



f'(R)[~g 2 r-U + H^ 



2H + 3H 2 



A-P A 



M 2 M 2 2M 2 



l + f(R))lg 2 ^+U + HiP) + 



+ 6(H + 3H 2 ) f'(R) g 2 ^- U + 



- Q[d d + 2Hd }f'(R) [g 2 ^ -i^ + H^ 
It is suitable to get, from system fl24|) -f)25l). the following equivalent one: 



R 



M 2 



4A - 3R$ + 9$ + 27m 



H= - 



2M 2 u 



where 



2(1 + f(R)) (g 2 ip 4 + + Htp) 2 ) - em -3$ + 3m 

= f(R) (g 2 ^ - (j, + Hip) 



r) 



(24) 



(25) 



(26) 
(27) 
(28) 



We can see that the term (1 + f(R)) [g 2 ^ 4 + (ip + Hip) 2 j corresponds to radiation since, if we 
neglect other terms, we get p = 3P. This result is a trivial generalization of the corresponding 
one in the model with minimal coupling (f(R) = 0), considered in |41j. In the f(R) modified 
model, T^ MS> has also terms proportional to f'(R). In this paper, we will show that these terms 
can actually play the role of the cosmo logical constant. 



3.3 Solutions with constant Hubble parameter 

We now investigate the de Sitter solutions for the model f[T9l . Our goal is to see how the 
YM field, which is described by £ym, can change the value of the cosmo logical constant. In 
particular, we will demonstrate in this section, that there do exist de Sitter solutions in the case 
when A = 0. 

We seek solutions with H = Hq =const, in other words, de Sitter and Minkowski solutions. 
If H = H , then R = R = 12 Hq, and (126 j) is a linear differential equation in 

+ 3H tf - AHl$ = B, (29) 

where the constant B = (MpR — 4A)/3. Eq. ( 12~9~|) has the following general solution 

i? = C x e m + C 2 e~ im - JL (30) 

and, from (!27j) . we get 

2(1 + f(Ro)) (VV 4 + (V> + ^) 2 ) - 3^ + = 0. (31) 
If fiYlH^) = — 1, then we have the equation 

$-H $ = 0, (32) 

which has the general solution: 

= C 3 + C 4 e Hot . (33) 
Thus, from (!30l) and (|33|) . we get that the de Sitter solution corresponds to 

B 

where C\ is an arbitrary constant. At A = 0, 

# dSo = C ie Hot - M 2 P . (35) 
It is easy to see that the Minkowski solutions (at j{12H'^) = — 1) correspond to 

$ M = C(t-t ), (36) 

where C and to are arbitrary constants. 

At = 0, Eqs. ^2M and (|25l) have the following nontrivial (i/; is not a constant) de Sitter and 
Minkowski solutions: 



^ = C ie **-^, (31) 
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• H = 0, A = 0, ip{t) = , /(0) = — 1, /x(0) is an arbitrary number. 

. H = H ? 0, A = 3M 2 P H 2 , j,(t) = ±§+Hoex % o{t _ to)) , /(0) = - 1, f (0) is an arbitrary 
number. 

These solutions do not change the value of the cosmo logical constant. Solutions, which corre- 
sponds to i9(t) ^ are more interesting. The following equation for ip arises 

f(R ) (f - s V + 2HW + /*V) = t|i " <V", (37) 
which, for C\ — 0, yields the following first order differential equation for ip: 

^ +2H ^=^'- H ^ 2 -iH!mwy (38) 

The trivial cases c = and H = have been considered above. In the general case, Eq. (|38|) 
does not satisfy the Fuchs conditions and, therefore, its solutions are multivalued functions (see, 
for example |44j). For C\ = the solution ip(t) can been found by quadratures, namely 



/ 1 " ■ ■ 



d ^ = t-t , c= ATT2finf)rT2y (39) 



~ g 2 r _ c ±H%f'(l2Ht) 
For nonzero values of C\ the solution ip(t) can be found numerically. 



4 Conclusion 

In this paper, we have studied cosmological models with an SU(2) Yang-Mills field. The model 
with a minimally coupled Yang-Mills field, described by the action ([[]), includes second- and 
fourth-order terms of the Yang-Mills field strength tensor. The second-order term can play the 
role of radiation, whereas the fourth-order one plays the role of the cosmological constant. 

We have shown that the function F{Z) can be reconstructed provided the Hubble parameter 
is given. In particular, we have demonstrate that de Sitter solutions exist only in the trivial 
case, namely when T{Z) is a constant. 

In order to obtain genuine de Sitter solutions, we have considered models in which the Yang- 
Mills field has a nonminimal coupling with gravity. We have explicitly shown that this model, 
described by the action in ( fl9|) . has de Sitter solutions even in the absence of a cosmological 
constant term. The de Sitter solutions correspond to the Yang-Mills fields which satisfy Eq. (13"5]) . 
This equation includes an arbitrary parameter. Depending on the value of this parameter, it has 
been shown that it can be easily solved in quadratures or, in the most general case, numerically. 
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Appendix 

In the FLRW spatially flat space-time 



T% = Ha 2 5 l3 , T) = H8), I* = 0, T% = 0, 

□ = - d d - 3Hd + (d 1 d 1 + d 2 d 2 + d 3 d 3 ) 

a{t) 2 



(40) 



#00 = - 3 ( H + H 2 j , R 0l = , Rij =[H + 3H 2 j 9ij , R = 6 [H + 2H 2 j . (41) 
Using the field A b in the form ([81) , we get two independent equations. The first one reads: 



-Roo — -Rgoo — 3H 



Ml vT °° 



(YM) 



+ A . 



(42) 



For any twice-differentiable functions f(R) and W(£), we get 

□ [f'(R)W] = -f"'(R)R 2 W - 2f"(R)RW - f'(R)W - SH f"(R)RW + f'(R)W 



g 00 a [f'(R)W\ - V do [f'(R)W\ = 3H f"(R)RW + f'(R)W 



Using 



1 



3 / 4? g 2 4- 4 



we obtain that Eq. ( 1421) is equivalent to 



2M 2 H 2 = (1 + f(R)) l^ + ^)-6(H + H 2 ) f'(R) < ' r ^ ' 



a? ' a 4 



2A 



a 2 / + 3 + 



+ 6H[Rf"(R)[^-^)+2f(R) l2 ~ 9 



7.2 A3( 



H4>) 



H<f>) 



(43) 



The second equation reads: 



Ru — -^Rga = — 9u[ + 3H 2 



Ml v n 



T (YM) + T (4) 



(44) 
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To calculate T^ YM ^ we use the following formulae (no summation over i) 



Fn F b n al3 — F b F b n 00 4- F b F b n jj 



1 _ 1 (g 2 ^ 



g^F^Ft-'-^F - 



2 V a 2 



+ 



and get (144p in the following form 

1 



2H - 3H Z 



2M 2 



1 + /(«)) f^ + ^!-2A + 



V a 



+ Q(H + 3H 2 )f'(R) ( ^ - £ ) - 6[^ + 2tf«9 ] (f{R) (^-^ 
Using (1T2]) . we rewrite Eqs. flS]) and (@ZD in terms of ip(t), to get ([24]) and (1251) . 



(45) 
(46) 



(47) 
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